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Abstract 

We have performed a joint analysis of prompt emission from four bright short 
gamma-ray bursts (GRBs) with the Suzaku-WAM and the Konus-Wind experiments. 
This joint analysis allows us to investigate the spectral properties of short- duration 
bursts over a wider energy band with a higher accuracy. We find that these bursts 
have a high E pea k, around 1 MeV and have a harder power-law component than that 
of long GRBs. However, we can not determine whether these spectra follow the cut-off 
power-law model or the Band model. 

We also investigated the spectral lag, hardness ratio, inferred isotropic radiation 
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energy and existence of a soft emission hump, in order to classify them into short 
or long GRBs using several criteria, in addition to the burst duration. We find that 
all criteria, except for the existence of the soft hump, support the fact that our four 
GRB samples are correctly classified as belonging to the short class. In addition, our 
broad-band analysis revealed that there is no evidence of GRBs with a very large 
hardness ratio, as seen in the BATSE short GRB sample, and that the spectral lag 
of our four short GRBs is consistent with zero, even in the MeV energy band, unlike 
long GRBs. Although our short GRB samples are still limited, these results suggest 
that the spectral hardness of short GRBs might not differ significantly from that of 
long GRBs, and also that the spectral lag at high energies could be a strong criterion 
for burst classification. 

1. Introduction 

The bimodal distribution of the duration of gamma-ray bursts (GRBs) indicates that 
there are two distinct classes of events. The long-duration bursts have typical durations of 
around 20 s, while the short- duration bursts (about 1/4 of the total) have durations of around 
0.3 s (Mazets et al. 1981; Norris et al. 1984; Dezalay et al. 1992; Hurley 1992; Kouveliotou et 
al. 1993; Norris et al. 2000). This suggests that the short and long GRBs may have different 
progenitors. Core collapse of massive stars and mergers of compact binaries are considered likely 
models for the long and short GRBs, respectively (Katz & Canel 1996; Ruffert Sz Janka 1999). 
Recently, thanks to the rapid position information provided by HETE-2 (Ricker et al. 2003) 
and Swift (Gehrels et al. 2004), afterglow observations have progressed dramatically and X-ray 
and optical afterglow emissions were discovered from some short GRBs (Hjorth et al. 2005; Fox 
et al. 2005), as well as long ones. Some short GRB afterglows have been found to be associated 
with galaxies not undergoing star formation, while some long GRBs have been found to be 
associated with energetic supernovae. These results support the hypothesis that long and short 
GRBs indeed have different progenitors. 

From the point of view of prompt gamma-ray emission, although the spectral charac- 
teristics of long GRBs have been well studied (Frontera et al. 2000: BeppoSAX; Kaneko et 
al. 2006: BATSE; Sakamoto et al. 2005b: HETE-2), our understanding of that of short GRBs 
is still incomplete, in part due to their very short durations. The spectral characteristics of 
BATSE short GRBs are often characterized by the hardness ratio, that is, the ratio of 100-300 
keV to 25-100 keV counts or fluence (Kouveliotou et al. 1993; Cline et al. 1999). These studies 
suggest that, although there is considerable overlap, short GRBs tend to be harder than long 
GRBs. Paciesas et al. (2001) and Ghirlanda et al. (2004) compared the spectral parameters 
of bright BATSE short GRBs with those of long GRBs by spectral fitting. They pointed out 
that the spectra of short GRBs are well described by a cut-off power-law model and that the 
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Band model (Band et al. 1993) did not improve the fit. They also confirmed that the spectra 
of short GRBs are harder than those of long GRBs. This was found to be a consequence of a 
flat low energy photon index, rather than a difference in the peak energies. 

There is another complicating issue in the classification of short and long GRBs. Since 
their duration distributions overlap, other distinguishing criteria have been proposed, such as 
differences in the host galaxy (Hjorth et al. 2003; Hjorth et al. 2005), spectral hardness (Cline et 
al. 1999), spectral lag (Norris 2002; Norris & Bonnell 2006), isotropic radiation energy (Amati 
et al. 2002; Amati 2006), and the existence of a soft hump (Norris & Bonnell 2006), or a 
combination of these and other characteristics (Donaghy et al. 2006). Indeed, some GRBs 
cannot be classified as short or long using the burst duration alone. For example, GRB 040429 
was classified as short using the burst duration (2.39 s), but other properties such as spectral 
lag and hardness resembled those of long GRBs (Wiersema et al. 2004; Fox & Moon 2004; 
Donaghy et al. 2006). GRB 051227 (Hullinger et al. 2005; Barthelmy et al. 2005a; Sakamoto 
et al. 2005) had a long enough duration of 8.0 s to be classified as long. However, this burst 
had no significant spectral lag in the initial spike, which is seen in many long GRBs (Norris 
et al. 1996), and it also had a long soft hump in the light curve. These two properties are 
similar to those of short GRBs. GRB 060614 (Gehrels et al. 2006) exhibited a long duration 
(102 s). However, it did not have any significant spectral lag in the initial spike, and also there 
were no signs of any associated supernova, despite its distance of z = 0.125, which should have 
been close enough to detect it. These examples show why the burst duration should no longer 
be considered a sole indicator for distinguishing between the short and long classes of bursts. 
Indeed Donaghy et al. (2006) have suggested the terms "short population bursts" and "long 
population bursts" to distinguish the two classes. In this paper, we will use simply "short" and 
"long" to describe these classes. 

Here, we report on the joint spectral analysis of four bright short GRBs simultaneously 
observed by the Suzaku-WAM (Suzaku Wide-band All-sky Monitor) and Konus-Wind. The 
WAM and Konus have wide energy ranges, from 50-5000 keV and 10-10000 keV, respectively, 
and the WAM has the largest effective area from 300 keV to 5000 keV of any current experiment 
with spectral capabilities. Thus, we can investigate the spectral characteristics of short GRBs 
up to MeV energies. We also classify the four GRBs as short or long using our data with high 
statistics up to the MeV energy region using a number of criteria: spectral lag, hardness ratio, 
isotropic radiation energy and existence of a soft hump. 
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2. Instruments and Observations 



2.1. Instruments 
2.1.1. Suzaku-WAM 

The Suzaku Wide-band All-sky Monitor (WAM) is the active shield of the Hard X-ray 
detector (HXD-II) (Takahashi et al. 2007; Kokubun et al. 2007) aboard Suzaku (Mitsuda et 
al. 2007). It consists of large-area, thick BGO crystals, and is also designed to monitor the entire 
sky from 50 keV to 5 MeV with a large effective area. The large effective area from 300 keV 
to 5 MeV (400 cm 2 even at 1 MeV) surpasses those of other currently operating experiments 
with spectral capability, and enables us to perform wide-band spectroscopy of GRBs with high 
sensitivity (Yamaoka et al. 2005; Yamaoka et al. 2006). The WAM is subdivided into four 
detectors, numbered WAM0 to WAM3, located at each side of the main detector of the HXD- 
II. We utilize the azimuthal (0 [degree]) and zenith (6 [degree]) angle to determine the incident 
direction of GRBs. The azimuthal zero angle ((f) = 0) is defined to be perpendicular to the 
WAM1 detector plane, and it increases in the order WAM0, WAM3, and WAM 2. The WAM0 
detector faces the solar paddle of the satellite and it always views the Solar direction. The polar 
axis that defines the zenith angle 9 is toward the Suzaku field-of-view of the main detectors. 
Thus the on-axis direction of each WAM detector correspond to (6, <fi) = (90, 0) for WAM1, 
(90, 90) for WAM0, (90, 180) for WAM3, and (90, 270) for WAM2. 

The WAM outputs two data types, the transient (TRN) data and the gamma-ray burst 
(GRB) data. The TRN data are always accumulated with 1 s time resolution and 55 energy 
channels. This can be used to monitor the bright soft gamma-ray sources with the Earth 
occultation method, as was done by CGRO/BATSE (Ling et al. 2000). On the other hand, the 
GRB data are recorded only for 64 s when the GRB trigger is activated, and the data cover 8 
s before and 56 s after the trigger time. The GRB data have four energy channels with 1/64 
s time resolution, in addition to the spectral data in 55 pulse height channels with 0.5 s time 
resolution. This allows us to perform both spectroscopy and timing analysis using the GRB 
data. In the initial phase before 2006 March 20, these two time resolutions for GRB data were 
set to 1/32 s and 1.0 s, respectively. 

The response matrix of the WAM is very complicated, depending on both the azimuthal 
and zenith angles of GRBs, because the WAM is inside the satellite body and suffers heavy 
absorption from it. In order to calculate the detector response of the WAM, we utilize a Monte 
Carlo simulation based on the Geant4 code, including many kinds of detector characteristics. 
In the pre-flight calibration, we measured these detector characteristics and developed a WAM 
response generator (Ohno et al. 2005). From the results of this pre-flight calibration, we found 
that the response generator can reproduce the measured spectral shape, but the measured 
effective area of the WAM varies drastically, depending on the gamma-ray incident angle. 
The trend of this angular response is different among the four WAM detectors. The angular 
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responses of the WAM1 and WAM3 detectors is complicated, and this trend is similar in both 
cases. We believe that this is caused by absorption from numerous electronics devices attached 
inside the satellite panels. The WAM2 detector is the most heavily affected, since the Ne/He 
chamber for the XRS is located inside the WAM2 field-of-view. This is not the case for the 
angular response of the WAMO detector, and it can be approximately described by a simple 
cosine function. This is because only the solar paddle is attached to the satellite panel, facing 
the WAMO detector. 

In order to estimate the uncertainty in the angular response, we performed an in-flight 
cross-calibration between the WAM, Konus, and the Swift/BAT, using the data of simultane- 
ously detected GRBs, and found that the measured spectral shapes agree quite well up to the 
MeV energy region, but that the observed flux fluctuates typically by 10-20% and by 40% at 
most above 100 keV. This uncertainty is consistent with pre-launch expectations. Below 100 
keV, there are still larger flux uncertainties of about 50%. We also confirmed that the WAMO 
detector has the most reliable detector response, with a flux uncertainty of 20% (Sakamoto et 
al. 2007, in preparation). 

2.1.2. Konus- Wind 

The Konus- Wind instrument (Aptekar et al. 1995) is a gamma-ray spectrometer aboard 
the GGS-Wind spacecraft. It consists of two identical detectors, SI and S2, which observe 
correspondingly the south and north ecliptic hemispheres in an all-sky monitoring mode. Each 
detector is a cylindrical Nal(Tl) crystal 13 cm in diameter and 7.5 cm in height. The experiment 
operates in triggered and waiting modes. In the triggered mode the burst time histories are 
recorded in three energy windows (Gl, G2, and G3) 0.512 s before and 229.632 s after the trigger, 
with a variable time resolution from 2 ms up to 256 ms. 64 energy spectra are measured in two 
partially overlapping energy ranges, nominally 10-750 keV and 0.2-10 MeV. In each range they 
are recorded in 63 channels with time resolutions starting at 64 ms, and subsequent spectra 
have resolutions from 256 ms to 8.192 s. In the waiting mode, each detector measures the 
count rate with a resolution of 2.944 s in three energy windows. For the GRBs considered here, 
the actual spectral ranges were 21 keV-16 MeV for the SI detector and 18 keV-14 MeV for 
the S2 detector; the actual time history windows were 21-83, 83-360, 360-1360 keV (SI), and 
18-70, 70-300, 300-1160 keV (S2). The detector response matrices depend only on the zenith 
(0) angle because of its axi-symmetric structure. 

2.1.3. The Interplanetary Network 

The Interplanetary Network (IPN) is a group of spacecraft with gamma-ray burst detec- 
tors. Its main purpose is to monitor the entire sky for gamma-ray bursts with close to a 100% 
duty cycle, and to localize bursts with up to arcminute accuracy. This is done by comparing 
the arrival times of a GRB among the spacecraft. Since the third IPN began operations in 
1990, over 25 spacecraft have participated in it, including Konus- Wind (Hurley et al. 2006a). 
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The Suzaku-WAM is the latest experiment to join the IPN. It now plays an important role 
in it, and many GRBs have already been localized using WAM data (Golenetskii et al. 2006a; 
Cummings et al. 2006; Hurley et al. 2006b; Golenetskii et al. 2006b). 

When a burst is detected by three widely separated spacecraft, a very accurate local- 
ization can usually be obtained. For the bursts discussed here, however, the Ulysses GRB 
experiment was off, and Mars Odyssey was the only interplanetary spacecraft. The resulting 
error boxes in this case tend to be large. Table 1 shows which spacecraft observed the events 
discussed in this paper. Note that GRB 060317 was not observed by Mars Odyssey, and there- 
fore the uncertainty in its localization is the largest. Table 2 gives the final coordinates of the 
error boxes. In most cases, the error boxes were not small enough and/or determined rapidly 
enough to warrant follow-up searches. In the case of GRB 060429, however, GCN Circulars 
were issued (Hurley et al. 2006b; Golenetskii et al. 2006b), and a Swift target of opportunity 
observation was carried out. However, for this and the other events, no counterparts were 
identified at any wavelength. 

2.2. The Short GRB Sample 

More than 70 GRBs triggered the WAM and Konus simultaneously between August 
2005 to October 2006. In order to perform joint spectral analysis of short GRBs, we selected 
well localized short GRBs and constructed detector response matrices. Only 8 short GRBs were 
localized well enough by Swift and/or the IPN in this period, and 4 of them (GRB 051127, GRB 
060317, GRB 060429, and GRB 060610) had hard enough spectra to detect in the highest energy 
band light curves of both instruments. Therefore, we select them to investigate the spectral 
properties of short GRBs up to the MeV energy region. 

GRB 051127 triggered the WAM and Konus at T (WAM) = 22:55:19.896 UT, and 
T (KW) = 22:55:15.860 UT 2005 November 27. The incidence angles on the detectors could be 
determined thanks to the IPN localization (Hurley 2005). The incident direction of this GRB 
was the nearest to on-axis for the WAM0 detector of the four, and thus the uncertainty in the 
detector response is the lowest for the WAM. GRB 060317 was very bright and it triggered the 
WAM and Konus at T (WAM) = 11:17:39.104 UT, and T (KW) = 11:17:35.996 UT 2006 March 
17. The incidence angle on the WAM was also good for the WAM0 detector. At T (WAM) = 
12:19:51.031 UT, and T (KW) = 12:19:49.712 UT 2006 April 29, the WAM and Konus detected 
the bright short GRB 060429. GRB 060610 triggered the WAM1 detector and the S2 detector 
of Konus at T (WAM) = 11:22:23.544 UT and T (KW) = 11:22:22.632 UT 2006 June 10. For 
this GRB, the incident direction was closer to on-axis for the WAM1 detector than WAM0. 
But this angle was determined to be good in the cross-calibration, and the flux uncertainty still 
should be within 20%. The detections are summarized in table 2. 
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3. Data Analysis 



We performed spectral analysis by applying three spectral models. The first one is a 
simple power-law (PL) model, 

where A is the normalization constant at 100 keV in photons cm -2 s _1 keV _1 , and a is the 
power-law photon index. 

The second model is a power-law with an exponential cutoff (CPL) model, 

where E pea k is the peak energy in the uF u spectrum, and represents the energy at which most 
of the power is emitted. 

The third model is a smoothly connected broken power-law model known as the Band 
model (Band et al. 1993): 

v ' V 100 / V V.k J (2 + a) 

VIOOy [ 100(2 + a) _ l Kf h ~ (2 + a) ' K J 

where a is the power-law photon index in the lower energy band, and (3 is the index in the 
higher energy band. 

First, we performed fits for each instrument alone. After confirming the consistency be- 
tween the WAM and Konus, we performed a joint fit. Selecting the time regions was important 
for this analysis because the WAM and Konus have different time resolutions for the spectral 
data. As mentioned in section 2, the WAM has a 1.0 or 0.5 s time resolution for spectral data, 
and can accumulate spectra both before and after the GRB trigger. On the other hand, Konus 
can accumulate spectra only after the trigger, and the time resolution varies depending on the 
burst intensity. In the case of the four short GRBs analyzed here, the time resolution of the 
Konus spectra is 64 ms for first four intervals and 8.192 s for remainder. Therefore, the time 
regions cannot be completely identical for the WAM and Konus, and we have selected the times 
for each instrument that are as close as possible between the two. To account for the differences 
in the accumulation times, we introduced a constant factor in the joint fits. The background 
spectra were extracted from both before and after the WAM time regions, but only after the 
time region for Konus. As mentioned in subsubsection 2.1.1, the detector response of the WAM 
strongly depends on the position of the GRB, and thus the error box of the IPN localization 
might cause some uncertainties for spectral analysis of the WAM. In order to estimate this 
effect, we calculated the detector response of the WAM for all corners of the box, as shown in 
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table 2, and performed spectral analysis by using each of them. This uncertainty is included 
as the systematic error. We used XSPEC version 11.3.2 for spectral analysis (Arnaud 1996). 
Throughout this paper, the quoted uncertainties are given at the 90% confidence level for one 
interesting parameter. A summary of the fitting results is given in table 3. 

We also examined the spectral lag using the cross-correlation function (CCF) between 
the light curves in two energy bands (Norris 2002; Norris & Bonnell 2006; Link et al. 1993; 
Band 1997; Yi et al. 2006). If fi(t) and V2{t) are the light curves in two different energy bands, 
the CCF is defined by 

CCF{T lag ) = ■ , where <r v = —\ i v i , (4) 

where N is the number of data points in the light curve. We used four energy bands for the WAM 
light curves: 50-110 keV, 110-240 keV, 240-520 keV, and 520-5000 keV. After calculating the 
CCF as a function of lag Ti ag , we obtained the peak value of ri ag by fitting it with a Gaussian 
profile (Yi et al. 2006). In this formula, positive spectral lag means that the spectrum has hard 
to soft evolution. 

3.1. GRB 051127 

The WAM and Konus light curves of GRB 051127 are shown in figure 1. The WAM 
light curve contains multiple spikes, and the T90 duration (the time to accumulate between 5 
and 95% of the counts), which is measured between 50 and 5000 keV, is 0.66 s. Thus, this 
burst is classified as short, using the burst duration alone. We extracted the spectrum from 
T (WAM) to T (WAM)+1.0 s for WAM and T (KW) to T (KW)+0.256 s for Konus, where 
the T s are summarized in table 1 and shown in figure 1. The propagation delay from Suzaku 
to Wind is 4.364 s for this GRB, i.e., correcting for this factor, one sees that the To(KW) = 
T (WAM)+0.328 s. We cannot adjust the time region further because the next Konus time bin 
has a duration of 8.192 s. The emission is clearly seen up to 5 MeV for the WAM and up to 2 
MeV for Konus, as shown in figure 2. We thus performed a fit from 100 keV to 5 MeV for the 
WAM and from 20 keV to 2 MeV for Konus. The WAM spectrum is well described by the CPL 
model with a photon index of — 0. 50^18 an d a peak energy E peak of 1257+^ keV. For Konus, 
the same model with a photon index of — 0.29to!2i an d an E pea k of 940 IlHq keV provides the 
best fit. The WAM spectrum seems to be slightly softer than that of Konus. We believe that 
this is because it includes a soft emission peak around T = To (WAM)— 0.2 s. There are no 
significant differences in \ 2 between the CPL and the Band model fits for either spectrum. The 
power-law photon index, a, in the low-energy portion and the peak energy, E pea k, are consistent 
with those obtained by the CPL model. However, we cannot constrain the power-law photon 
index, /3, at high energies and obtain only an upper limit of /3 < —2.03. We then performed a 
joint fit. The CPL model with photon index — 0.44lo;]J and E pea k 1168^i5o keV provides the 
best fit. Table 3 summarizes the results of all of the conducted fits. 
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3.2. GRB 060317 



There are two intense pulses in the light curve of GRB 060317, as seen in figure 1. 
The burst has a total duration of about 1.31 s for the WAM. We cannot distinguish whether 
this GRB belongs to the short or long population from the burst duration alone, because the 
duration distribution overlaps around 1-2 s. In order to extract the spectrum from both pulses, 
the WAM time region was taken to be T (WAM) to T (WAM) +4.0 s. For Konus, we had to 
select a larger time region than that of the WAM, so as to include both pulses: from T (KW) 
to T (KW)+8.336 s. Konus detected this burst 3.49 s before the detection by the WAM, and 
thus T (KW) corresponds to T (WAM)+0.389 s. We can see that emission is clearly detected 
up to 5 MeV for WAM and up to 10 MeV for Konus, as shown in figure 2. Therefore, we 
performed the spectral fitting from 100 keV to 5 MeV for the WAM and from 20 keV to 10 
MeV for Konus. We cannot constrain E pea k from the WAM spectrum, and a simple power-law 
with a photon index of — Lillys fits the observed spectrum well. On the other hand, the 
Konus spectrum is well described by the CPL model with a photon index of — 1.02lg!i3 and an 
Epeak of 2089tg34 9 keV. The Band model does not provide any improvement to the fit. When 
we performed the joint fit, E pea k became very large (5687i?jg§ keV). This might come from the 
large IPN error box, which is about 5 deg 2 , and the systematic uncertainty might be larger than 
our estimation. Therefore, we take the spectral parameters of this burst to be those obtained 
by Konus. 

3.3. GRB 060429 

This GRB shows a very simple short spike structure with a duration of 0.08 s in the 
WAM light curve (figure 1); based on duration alone, it is short enough to be placed in the 
short class. Since we changed the WAM time resolution from 1.0 s to 0.5 s for spectral data, 
and from 1/32 s to 1/64 s for light curve data on 2006 March 20, we could select the time 
region to be T (WAM) to T (WAM)+0.5 s, and obtain a better S/N. For Konus, we selected 
the time interval from T (KW) to T (KW)+0.128 s. The propagation delay from Suzaku to 
Wind is 1.530 s for this GRB, i.e., correcting for this factor, one sees that the To(KW) = 
T (WAM)+0.211 s. We can see in figure 2 that the emission extends up to 5 MeV for both 
WAM and Konus. We thus performed spectral fittings from 100 keV to 5 MeV for the WAM 
and from 20 keV to 5 MeV for Konus, respectively. These spectra are also well described by the 
CPL model with a high E peak . The photon index is -0.99tg;iJ for the WAM and — 0.76±g;i| for 
Konus, respectively. The E pea k can be constrained to 1583137? keV h J Konus, but the WAM 
data do not constrain it. The fit does not improve in the Band model. The joint fit also gives 
consistent results. 
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3.4- GRB 060610 

We can regard this GRB as short, based on the WAM burst duration of 0.79 s alone. 
The light curve exhibits a double-peak structure, as shown in figure 1. The time regions for 
spectral analysis were selected from T (WAM)+0.5 s to T (WAM) + 1.0 s and from T (KW) 
to To(KW)+0.256 s. The propagation delay between these two instruments was 1.402 s, i.e., 
correcting for this factor, one sees that the To(KW) = To (WAM) +0.490 s. Figure 2 shows 
that the emission extends up to 5 MeV for both instruments. We thus performed spectral 
fittings from 100 keV (WAM) and from 20 keV (Konus) to 5 MeV. Both spectra are also 
well described by the CPL model. The WAM and Konus photon indices are — 0.90ig; 2 2 and 
-0.94lo;ii, respectively. The E peak is constrained by the WAM and Konus to 28891^75 keV 
and 1536l|gg keV, respectively. From the joint fit, we obtained a photon index of — 0.92i ; 12 
and an E peak of 1821±|§t keV. There is no significant improvement of the fit with the Band 
model, and we can not constrain the power-law photon index, (3, at high energies; we obtained 
an upper limit of j3 < —2.50. 

3. 5. Spectral Lag 

Figure 3 shows an example of the cross correlation function (CCF) calculated for GRB 
060429. The CCFs of our four short GRBs all exhibit a sharp shape structure around zero. 
We estimate the peak values of the CCFs using Gaussian profiles (Yi et al. 2006). The results 
are given in table 6. The spectral lags are all consistent with zero for all energy ranges. For a 
comparison, we also calculated the CCF of three bright long WAM GRBs: 050924, 060915, and 
060922. Their CCFs display broad peaks, as shown in figure 3. Thus, we used a polynomial 
instead of a Gaussian profile to estimate the peak values; the results are also given in table 
6. We found that these long GRBs have significant positive lags, which become larger towards 
the higher energy bands. We confirmed that this lag is not caused artificially, by analyzing 
simulated data with no spectral lag, and indeed obtained zero lag. 

4. Discussion 

4-1- Spectral Characteristics of Short GRBs 

The spectra of all four GRBs are well described by a power law with an exponential 
cut-off (CPL), and the fits do not improve significantly when we apply the Band model. To 
investigate the possibility that the spectra are better described by the Band model, we used 
the Band model with a typical value of the power law photon index in the high-energy portion 
of (3 — —2.3. When we performed joint fits with this model, the x 2 /d.o.f became 66/55, 69/58, 
and 107/75 for GRB 051127, 060429, and 060610, respectively. For GRB 060317, there was 
also no improvement in the Konus spectrum with this model, whose % 2 /d.o.f is 53/49 (see table 
3). Therefore, we cannot determine whether these GRBs follow the CPL model or the Band 
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model from our data. The power-law photon index, a, of the low-energy part of our sample is 
distributed around —1.0 and the mean value is —0.73. This is larger than that of the BATSE 
long duration GRBs of —1.05 (Ghirlanda et al. 2002). In particular for GRB 051127, there 
is a significantly flatter photon index of —0.41, which is larger than —2/3. This distribution 
of a is also reported for BATSE bright, short duration GRBs (Ghirlanda et al. 2004). They 
showed that the spectra of their short GRB sample were also described by a CPL model. If a 
flat slope is a global characteristic of the prompt emission spectrum of short GRBs, the simple 
synchrotron emission model may no longer be applied and other emission mechanisms, such as 
jitter radiation (Medvedev 2000), are required. 

4-2. Classification of Short GRBs 

It is ambiguous in some cases to classify GRBs using the burst duration alone, because 
the observed long and short duration distributions overlap around 1.0-2.0 s. The dilation of 
the intrinsic duration of GRBs due to cosmological distance, and the dependence of duration 
on energy, must also be considered. Therefore, we utilize several other methods of classifying 
short and long GRBs below. In this section, we discuss whether the four bursts are really in 
the short class by examining these criteria: 1) spectral lag, 2) spectral hardness, 3) radiation 
energy, and 4) extended emission. So far, these criteria have been confirmed mainly based on 
the BATSE data. We can test these criteria in a higher energy band than BATSE by the WAM 
and Konus data with high statistics even in MeV energy band. 

The properties of the spectral lag are reported to differ between short and long-duration 
GRBs. Short GRBs exhibit lags, which are consistent with zero, while long GRBs have sig- 
nificantly non-zero positive spectral lag (Norris et al. 2000; Norris & Bonnell 2006). We have 
measured the spectral lags using the WAM light curves. Figure 4 plots the relation between 
the duration and the spectral lag for four short GRBs, for 50-110 keV vs. 110-240 keV, 50-110 
keV vs. 240-520 keV, and 50-110 keV vs. 520-5000 keV. We also plot the lags of three bright 
long duration WAM GRBs for comparison. We find that the spectral lag of all four short 
GRBs is consistent with zero in approximately the same energy bands as BATSE (50-110 keV 
vs. 110-240 keV or 50-110 keV vs. 240-520 keV), while non-zero lag, evolving with T90, is 
seen for the three long GRBs. Therefore, the WAM data confirm the difference in the spectral 
lag obtained with BATSE. In addition, figure 4 demonstrates the WAM capabilities in the 
high-energy band: lag evolution is seen above 520 keV for long GRBs. 

The hardness ratio of short GRBs is reported to be slightly harder than that of long 
GRBs from BATSE results (Cline et al. 1999; Paciesas et al. 1999). Although the boundary 
between short and long GRBs in the duration-hardness plane is ambiguous, we nevertheless 
consider this difference as a possible criterion. Therefore, we obtained the spectral hardness 
of our sample of short GRBs, based on the best-fit spectral model. We derived the hardness 
ratio using the 50-100 and 100-300 keV fluence, which is same energy range as BATSE. We 
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plot them against the T90 duration in figure 5, together with the BATSE results (Paciesas et 
al. 1999) The hardness ratios of our short GRBs are distributed from 3.61 to 5.64 with an 
average value of 4.61, and all have greater hardness ratios than the average value of the BATSE 
long (T90 > 2.0) GRBs (3.27). Therefore, all four are consistent with a short classification. 
However, these short GRBs do not show a clear separation from the long GRBs, as in the 
BATSE samples in figure 5. Although all our short GRBs have an E pea k around 1 MeV, none of 
them exceeds a hardness ratio of 8, whereas half of the BATSE short GRBs do. This is because 
the hardness ratio depends on the low-energy photon index rather than E pea k (Ghirlanda et 
al. 2004; Sakamoto et al. 2006). The short GRB which has the highest hardness ratio in our 
sample is GRB 051127, with 5.64 tl'lt- GRB 051127 also has the hardest photon index, -0.44, 
but its Ep ea k is one of the lowest in our sample. As mentioned in Sakamoto et al. 2006, the 
low-energy photon index has to be flatter than to reproduce a hardness ratio greater than 8, 
which is seen in the BATSE "very hard" short GRBs. Based on our broad band spectroscopy 
of bright, high E pca k short GRBs, none has a hardness ratio similar to the BATSE "very hard" 
short GRBs. 

For long GRBs, a correlation between the peak energy and the isotropic equivalent 
radiation energy in the source frame has been reported: the Amati relation (Amati et al. 2002). 
However, short GRBs whose redshifts are known have lower isotropic energy than long GRBs 
and do not satisfy this relation (Amati 2006). We calculated the isotropic energy, E iso , of our 
four short GRBs from 1.0 to 10000 keV in the source frame, assuming redshifts from 0.1 to 
10.0. We used cosmological parameters (H , fl\, fl m ) = (65, 0.7, 0.3). We included the effects 
of redshift on the energy spectra and E peak of these bursts. Figure 6 shows the results of the 
calculation in the E p -Ei SO plane together with the HETE-2 and BeppoSAX results (Sakamoto 
et al. 2005b). We could not find any redshift that satisfies the Amati relation. All of the bursts 
in our sample exhibit a lower isotropic radiation energy by more than one order of magnitude 
when compared with long bursts of the same peak energy. This is consistent with previous 
results. 

Long, soft emission humps a few tens of seconds after the short main spike were reported 
for several short GRBs (Norris & Bonnell 2006), particularly GRB 050709 (Villasenor et al. 
2005) and 050724 (Barthelmy et al. 2005a), which have measurements of the redshift and 
detections of the host galaxy. Thus the existence of the soft emission hump is considered to be 
one of the criteria to classify short GRBs. We searched for long soft humps in the time histories 
of our four short GRBs, but did not detect them. To do this, we calculated the upper limits 
to the peak soft photon flux at the 3 a level for the WAM and Konus. In this calculation, 
we assumed that the spectrum of the soft emission was a simple power law with photon index 
-2.0. We found 3 a upper limits to the 2.0 to 25 keV peak flux of 4.47, 6.63, 9.21, and 7.54 
photons s" 1 cm" 2 , for GRB 051127, 060317, 060429, and 060610, respectively, using both the 
WAM and Konus. All of these values exceed 2.72 (± 0.47), which is obtained for GRB 050709 
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(Villasenor et al. 2005). Therefore, the non-detection of the soft emission for these bursts adds 
no constraint: the extended emission could be too soft to be detected by the WAM and Konus. 

5. Conclusions 

We performed a joint analysis of four bright, hard-spectrum, short-duration GRBs, lo- 
calized by the IPN, with the Suzaku-WAM and Konus- Wind, in order to investigate the spectral 
properties of short GRBs into the MeV energy band. From the spectral analysis, we found that 
these bursts have a high E pea k around 1 MeV; the spectral parameters can be constrained tightly 
by joint fitting. However, we could not determine whether these spectra follow the CPL model 
or the Band model, because there is no significant improvement between these models. The 
power-law photon index is slightly flatter than that of long BATSE GRBs. In particular, GRB 
051127 has a very flat spectrum with a photon index, a, that is larger than the synchrotron 
limit of -2/3. 

We also examined several other spectral and temporal properties to confirm that the 
bursts belong to the short class. Most properties, such as the spectral lag, the spectral hardness, 
and the isotropic radiation energy, indicated that the four bursts were indeed in this class. In 
addition, a broad-band spectral analysis with the WAM and Konus revealed some interesting 
properties of short GRBs: i) There is no evidence for very hard GRBs that have the large 
hardness ratios seen in BATSE short GRBs. ii) The spectral lag of our short GRBs does 
not show any clear evolution, even in the MeV energy region, unlike that of long hard GRBs. 
Although our short GRB sample is still limited, these results might suggest one possibility, that 
the hardness ratios of the short GRBs are not so different from those of long GRBs, and also 
that the spectral lag at high energies can be a strong criterion for classifying short and long 
GRBs. 
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Fig. 1. Background-subtracted WAM (top) and Konus (bottom) light curves of our sample of short GRBs. 
The zero times of these light curves are adjusted to the Konus trigger time. The time interval for extracting 
the spectrum is shown by the dashed-line. Left-top panel is GRB 051127, right-top panel is GRB 060317, 



Table 1. Interplanetary Network Spacecraft Observations of the four WAM and Konus short GRBs. 



Trigger time (UT) 
Suzaku-WAM Konus-Wind 


Spacecraft*^ 


GRB 051127 
22:55:19.896 


22:55:15:860 


Konus, Suzaku, Odyssey HEND, HETE-FREGATE 


GRB 060317 
11:17:39.104 


11:17:35.996 


Konus, Suzaku, INTEGRAL SPI-ACS 


GRB 060429 
12:19:51.031 


12:19:49.712 


Konus, Suzaku, Odyssey HEND & GRS, INTEGRAL SPI-ACS, RHESSI 


GRB 060610 
11:22:23.544 


11:22:22.632 


Konus, Suzaku, Odyssey HEND & GRS, INTEGRAL SPI-ACS, RHESSI 



*: HEND and GRS arc the High Energy Neutron Detector and Gamma-ray Sensor head (Hurley et al. 2006a) 
f: SPI-ACS is the spectrometer anti-coincidence detector (Rau et al. 200.5) 
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Fig. 2. Joint spectral fitting with a cut-off power law (CPL) model for the WAM (square) and the Konus 
(triangle) data. The solid lines represent the best-fit CPL model. 
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Table 2. Error box centers and corners of IPN localization and Incident angle. 



GRB 


Area, sq. deg. 


o^ooo.o 


^2000.0 




Incident an§ 
6, <p (WAM) 


de [degree]* 
9 (Konus) 


051127 


0.13 


301.500 


-37.782 


(center) 


56.3, 80.2 


17.4 






301.371 


-36.873 




56.6, 79.1 








301.695 


-38.523 




56.4, 78.6 








301.622 


-38.663 




56.0, 81.2 








301.302 


-37.024 




56.6, 79.3 




060317 


5.3 


301.900 


-10.882 


(center) 


125.0, 49.5 


9.1 






301.391 


-17.546 




131.7, 49.5 








301.976 


-8.163 




122.3, 49.4 








303.343 


-4.329 




118.4, 50.7 








301.987 


-13.575 




127.7, 49.9 




060429 


0.14 


115.314 


-24.926 


(center) 


146.1, 134.1 


45.5 






115.307 


-24.932 




146.1, 134.1 








113.108 


-25.520 




145.5, 130.5 








117.509 


-24.233 




146.5, 137.8 








115.321 


-24.919 




146.0, 134.0 




060610 


0.013 


354.762 


52.107 


(center) 


41.1, 23.6 


48.3 






354.294 


52.375 




40.7, 23.5 








354.260 


52.417 




40.6, 23.5 








355.227 


51.835 




41.4, 23.7 








355.261 


51.793 




41.5, 23.7 





*: 6 and <fi are the zenith and azimuthal angles, respectively, in degrees. 
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Table 3. Results of spectral fittings. 



GRB 051127 



Model 


Detector 


a* 






Epoak [keV] 


Norm 1 '- 


Const § 


X 2 /d.o.f 


PL 


WAMO 
KW(S1) 


— 1 94+ - 05 

— 1 07 + 






- 


q 7O+0.41 
°- ' °-0.36 

4 fjy+o-o 8 

^■°'-0.59 




142/24 
70/31 




WAM0+KW(S1) 


1 - zz -0.04 








°- ou -0.54 


o.67l°; u7 


229/56 


CPL 


WAMO 


U - JU -0.18 






1257t?«° 


Z - JO -0.35 




29/23 




KW(S1) 


-0.291°!? 








r q 7 +0.87 
-°' -0.82 




32/30 




WAM0+KW(S1) 


-o.44±8-.ig 






nests 


5.181°;*? 




64/55 


Band 


WAMO 


-o-3o±g-iS 


< 


-2.08 


9991^ 


2.48±g:S 




32/22 




KW(S1) 


-0.34l°f 7 


< 


-2.03 


qqq+58-5 
yyy_i87 


r on+0.87 
O-OU_ gg 




32/29 




WAM0+KW(S1) 


-0.321°;^ 


< 


-2.33 


10381^ 


5 1fi +0 ' 58 
°- 1D -0.56 


o.48l o: ° i 


64/54 


GRB 060317 


Model 


Detector 


a 







Epcak 


Norm 


Const 


X 2 /d.o.f 


PL 


WAMO 
KW(S2) 


i 1-1+0.05 
i - 11 -0.06 

1 OO + 0.05 

1 - oz -0.Q5 








2 14+ - 48 

z - 14 -0.34 

Q0+ - 08 

u - yu -0.08 




54/25 
74/50 




WAM0+KW(S2) 


-1 1fi +005 

J— 10 -o.o5 








u - ' u -0.08 


q 97+0.25 
J - z ' -0.24 


167/76 


CPL 


WAMO 


1 qo+0.05 
1 - uo -0.06 






>5000 


2 06+ 48 

z - UD -0.33 




50/24 




KW(S2) 


-1 02+ 015 






20891^° 


n 07+0.09 

U - 8 '-0.09 




53/49 




WAM0+KW(S2) 


-i-03±8:g 






5687lf^ 


70+ - 08 

u - ' u -0.08 


q r)7+0-25 
J - u ' -0.22 


144/75 


Band 


WAMO 


-0.991°^ 


z - ou -0.09 


>5000 


no +0.48 
z - uz -0.33 




47/23 




KW(S2) 




< 


-2.73 


qqy+1089 
aa ' -1281 


OQ+0.08 
u - ao -0.10 




55/48 




WAM0+KW(S2) 


-o.99±8:2S 


< 


-1.58 


qnqn+6743 


76+ - 03 

u - ' u -0.14 


9 74+0.26 

z - '^-0.19 


128/74 



*:power law photon index for PL and CPL model, and low energy index for the Band model. 
f:powcr law photon index for the high energy component of the Band model, 
^normalization in photons keV _1 cm -2 s — 1 at 100 keV. 

§:constant factor for the WAM to Konus normalization in joint spectral fittings. 
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Table 4. (Continued.) 



GRB 060429 


Model 


Detector 


a 


P i^pcak 


INorm 


Const 


2 / j r 

X /d.o.j 


PL 


T T T A "1\ IT f\ 

WAMO 


-!-24 lo.io 




3.66+ 7 ° 




42/25 




T7~~\ T7/0 1 \ 

KW(bl) 


-i ot+0.04 

- 1 -27I . 4 




q rr+0.85 

9 - 00 -0.85 




98/33 




\\ T A A If Fl i ~\Z \ H T ( C 1 '\ 

WAMU+KW(bl J 


1 o«+0.05 




9 - 50 -0.85 


n ,m+0.05 


14U/59 


CPL 


WAMO 


-0.99 ll 2 l 


>1430 


, | .1/1 -1- 1 ^ 1 

°- za -0.83 




37/24 




KW(S1) 




1583t 4 ^ 


11 + 117 

-1.14 




30/32 




WAM0+KW(S1) 


-o.aolg : g 


1 720+ 504 

1 ' zu -384 


10 9+F 3 


o-28±8:S 


70/58 


Band 


WAMO 


_n qq+0-55 


< -1.28 >1000 


o 90+0-78 




37/23 




KW(S1) 


-0 62 + 14 

U - Dz -0.29 


< -2.02 1560t 4 ^ 


n o +1 - 25 

ii - u -1.09 




30/31 




WAM0+KW(S1) 


-0 62+ 12 

U - Dz -0.17 


< -1.67 1251±|g? 


i i o +0.97 
11 - z -1.10 


0.27t°; o ^ 


73/57 


GRB 060610 


Model 


Detector 


a 


P Epcak 


Norm 


Const 


X 2 /d.o.f 


PL 


WAM1 


-1 16+ 008 

1 - iu -0.09 




o 4 6 +0.82 
°-^ u -0.73 




38/25 




KW(S2) 


-, oq+0.04 
1 - oo -0.04 








123/50 




WAM1+KW(S2) 


-1 29+ - 04 

1 - ZJ -0.05 




5-98 j£g 




175/76 


CPL 


WAM1 


-o.9oig:g 


2889t?^g 


o i i +0.83 
J ' il -0.74 




24/24 




KW(S2) 


-0.94ig : g 


1536i^ 


7 in+0.66 




67/49 




WAM1+KW(S2) 


-0.92j£g 




6 71 + - 63 


u - oz -0.09 


104/75 


Band 


WAM1 


-0.81±°il 


< -8.82 1799l^f 


o i 9+0.83 
°- iz -0.73 




24/23 




KW(S2) 


-0.99i°;SI 


< -2.21 1531±l?f 


y on+0.58 
' - zy -0.73 




69/48 




WAM1+KW(S2) 


-0.88l°i° 


< -2.50 1687tg 4 


6.70t°;^ 


o 5n+ 016 

u - ou -0.07 


102/74 



Table 5. Fluence and peak flux of our four short GRBs. 



GRB 



Fluence* 



Peak fluxt 



[10 6 erg cm 2 ] 
50-100 keV 100-300 keV 300-500 keV 500-1000 keV 



[10 6 erg s 1 cm 2 ] 



[ph s 1 cm 



GRB 051127 
GRB 060317 
GRB 060429 
GRB 060610 



+0.02 
-0.02 

+0.06 
0.06 

+0.01 
0.01 



13+ 001 
u - io -0.01 



0.16 
0.72 
0.10 



0.93 



+0.08 
-0.09 



9 60+ 002 

z - DU -0.02 



0.49 



+0.04 
0.04 



^+ - 04 

U - OO -0.04 



+0.11 
-0.11 

+0.02 
0.02 

+0.05 
0.05 



^2+ - 05 

U ' OZ -0.05 



1.1 

2.38 
0.49 



2.31 
5.47 



+0.24 
-0.24 

+0.04 
0.04 



1 10+ 012 
1 - 1U -0.12 

1 12+ 011 

i ' lz -0.10 



28 



+3.2 
-3.2 



67±I;? 

qq+2.3 
°°-2.3 



22 



+ 2.5 

-2.5 



41+3-4 

^-3.3 



77 



+7.7 
-6.6 



04+2.3 
°^-2.3 



The energy fiuence obtained by the CPL model in the joint fitting of the WAM data. 



f: The energy and photon peak flux in 1/32 s for GRB 051127 and 060317, and in 1/64 s for GRB 060429 and 060610, 
respectively, obtained by the WAM 50-5000 keV data. 
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Table 6. Duration(Tgo) and spectral lag of four short GRBs and two long GRBs obtained by the WAM data. 



GRB Date 


T 90 * 7 

[s] 


•lag (TH0,TH1 


t T! ag (TH0,TH2)t 
[ms] 


r lag (TH0,TH3)t 








short GRBs 




051127 


0.66 


17 +18 
1 ' -21 


io + 14 
J-O-14 


7 i +23 
' • 1 -18 


060317 


1.31 






7 o+8.4 
' -°-7.2 


060429 


0.08 




12+ 14 

iz -15 


11+26 
11 -28 


060610 


0.79 


i 7+33 
11 -25 


oi +36 
zi -50 


9D+ 36 
zu -45 








long GRBs 




050924 


29.09 


2631" 




1 OKE +525 
lz00 -266 


060915 


43.16 


B6tg 




i7o; 7 ^ 


060922 


13.88 


199+41 


505±1« 


1004±ii? 



*: T90 duration measured in 50-5000 kcV energy band of the WAM light curve. 

f: The energy bands of TH0, TH1, TH2, and TH3 correspond to 50-110 kcV, 110-240 kcV, 240-520 keV, and 
520-5000 keV, respectively. 




-1 -0.5 0.5 1 -10 -5 5 10 

Lag (s) Lag (s) 

Fig. 3. Example of the light curves of the WAM data (top panel) and the CCFs (bottom panel) for 
a short GRB 060429 (left) and a long GRB 060922 (right). The solid line shows CCF(50-110 keV to 
110-240 keV), the dashed line shows CCF(50-110 kcV to 240-520 kcV), and the dot-dashed line shows 
CCF(50-110 keV to 520-5000 keV). 
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Fig. 4. Relation of the Tgo duration and the spectral lag. From top to bottom, the lag of 50-110 keV to 
110-240 keV, 50-110 keV to 240-520 keV, and 50-110 keV to 520-5000 keV arc shown. The filled circles 
are for the four short GRBs of this paper and the crosses are for long GRBs. 
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Fig. 5. Relation between the duration (Tgo) and the hardness ratio of 100-300keV to 50-100 keV. The 
Dots show BATSE results. The results of our four short GRBs are shown as filled squares. 




Fig. 6. Relation between E poa k and Ej so for our four short GRBs (filled square. GRB 051127; filled circle, 
GRB 0603f7; filled triangle, GRB 060429; and open diamond, GRB 0606f0). The trajectories were 
calculated by assuming redshifts from 0.1 to 10.0. The data points along the trajectories correspond to 
z=0.1, 0.2, 0.3, 0.5, 0.7, 1.0, 2.0, 3.0, 5.0, 7.0, and 10.0, respectively. The dasheddine represents the Amati 
relation E pca k, sr c oc EP^ for long GRBs. 
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